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Abstract 

We have calculated the Wandzura-Wilczek contribution to the twist-3 part of 
7*7 — ► 2-7T amplitude. It describes interaction of the longitudinally polarized virtual 
photon with the real one, and it is suppressed by 1/Q, where Q 2 is the virtuality 
of the 7*, as compared to the twist- 2 contribution. We have found that, in the 
Wandzura-Wilczek approximation, factorization applies to the twist-3 amplitude. 

1 Introduction 

Hadron production in the reaction j*(q)j(q') ~ > hadron(s), q 2 = —Q 2 , has been a sub- 
ject of considerable interest in the QCD community from both experimental |], § and 
theoretical [||, H |5j points of view. Recently it has been proposed [|, [7| to investigate a 
process 7*7 — > titt when the two pion state has a small invariant mass. Thanks to the 
QCD factorization theorem ||, in the leading twist approximation the amplitude can be 
represented as a convolution of perturbatively calculable Wilson coefficients and new non- 
perturbative objects, the so-called two-pion distribution amplitudes (27rDA's). They are 
given by matrix elements of twist- 2 QCD string operators between vacuum and the two- 
pion state 0, ||. Moreover, 27rDA's can be related by the crossing symmetry to skewed 
parton distributions [pi], |ll|] which recently have been subject of considerable interest. 



Later, NLO corrections to the leading twist amplitude have been calculated and the 
recent thorough phenomenological analysis |13| have shown that experimental studies of 
27r production cross-section are possible with existing e + e~ facilities. 
In general, parity invariance restricts the number of independend helicity amplitudes in 
the 7*7 — > 7T7T process to three, which correspond to different possible projections of the 
total angular momentum on the 7*7 collision axis, J z = 0, ±1, ±2. J z — and J z = ±2 
amplitudes scale as Q 2 — > 00, and their corresponding QCD representation has been 
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discussed in [0, |9], |T^]. The amplitude corresponding to J z = ±1 is suppressed as 1/Q and 
has not been discussed so far. On the other hand, as it has been demonstrated in it 
can be extracted from the experimental data by a suitable weighting procedure. 
The purpose of this paper is to fill this gap in the QCD description of the 7*7 — > irn am- 
plitude. As it can be anticipated from its 1/Q suppression, the leading contribution to the 
J z = ±1 amplitude is given by matrix elements of twist-3 quark-quark and quark-gluon 
operators. In this paper we restrict ourselves to the Wandzura-Wilczek (WW) approx- 



imation |TJ] i.e., we consider only the contribution from quark operators. Recall that 
the recent experimental analysis of the polarized nucleon structure function gr{x, Q 2 ) 



indicates that the WW approximation can account rather well for the experimental data. 
A similar conclusion is reached in the instanton model of the QCD vacuum |16), where 
quark-gluon correlations in a nucleon are suppressed by a small parameter given by the 
packing fraction of instantons. Our results can be therefore considered as an estimate of 
the order of magnitude of the J z = ±1 amplitude. 

The problem of twist-3 contributions to the hard exclusive amplitudes has recently ac- 
quired a considerable attention [[H], ITS], [H5|, BT], |22fl . In particular, the approach pre- 



sented here is rather similar in spirit to the analysis of pi] , which presented a detailed 
consideration of the problem of the DVCS amplitude on a pion. 

2 General definitions 

Kinematics of the reaction 7*(g , )7(<7 / ) — > 11^1)71^2) can conveniently be described in 
terms of a pair of light-like vectors p, z which obey 

p 2 = z 2 = 0) p . z ^ (1) 

and define longitudinal directions. Here p ■ z = p^ . Let P and k denote total and 
relative momenta of the ir meson pair, respectively, 

P2 = + = W 2^ k 2 = fa _ k ^2 = 4m 2^_ W 2^p. k = Q ( 2 ) 

The initial and final states momenta can be decomposed as 

Q 2 2 n 2 , Q 2 + w 2 /2 

2{p ■ z) 2{p ■ z) 

W 2 

P = q + q ' = p + Z) P 2 = w 2 (3) 

2(p ■ z) 

k = £p- — r z + k ± 

2(p ■ z) 

The longitudinal momentum distribution between pions is described by the variable £ = 
(k ■ z)/(p ■ z). Alternatively, 

£ = (3 cos 9 cm , 

where 8 cm is the polar angle of the pion momentum in the CM frame with respect to 
the direction of the total momentum P and (3 is the velocity of produced pions in the 
center-of-mass frame 



Ami 



The amplitude of hard photo-production of two pions is defined by the following matrix 
element between vacuum and two pions state: 



T"" = ij d 4 xe- ix ^(27c{P, k)\TJ"{x/2)J v {-x/2)\Q) , q=iz{q- q') (4) 

where J M (x) denotes quark electromagnetic current. Hard photo-production corresponds 
to the limit Q 2 S> W 2 > Aq CD where the amplitude (f|) can be represented as an expan- 
sion in terms of powers of 1/Q. According to the factorization theorem the leading twist 
term in the expansion can be written as a convolution of hard and soft blocks. The coef- 
ficient functions can be calculated from appropriate partonic subprocesses 7* + 7 — > q + q 
or 7* + 7 —>■ g + g. Soft blocks, describing transition of partons into the final meson pair 
are given by quark and gluon 2n distribution amplitudes. 



According to the analysis of Ref. [12] the amplitude T^ v is a sum of three terms 



T^{q, q', P, k) = l ~{-gn T TZ™{q, q' , P, k) + i^A^P + <J)"TF r {q, q' , P, k) + 

+ l^VTM,P,k) (5) 

where {—g^ u )T — { p '" z p. z z * — 9^ is the metric tensor in the transverse space and k^fk 1 ^ 
denotes traceless, symmetric tensor product of relative transverse momenta 
Amplitudes To" r7r and T^™ correspond to J z = and J z = ±2, respectively. They have 
been discussed at length in |7|, |8|, || [12] In this paper we consider the J z — ±1 amplitude 
T^ nn . It arises when the (real) transverse photon collides with the (virtual) longitudinal 
one and produces a pair of pions in L z = ±1 state. It is easy to see that such a process 
requires helicity flip along the quark line, and hence it vanishes in the leading twist 
approximation. 

Consider now a frame in which the pion pair is moving in the positive z direction so 
that p + and z~ are the only nonzero components of p and z, respectively. In an infinite 
momentum frame p + ~ Q — > 00 with fixed (p • z) ~ 1. From (|3|) it follows that in this 
frame (k ■ z) ~ 1 and k± ~ Q°. Hence, the contribution due to X^ 71 " 71 " is suppressed by 1/Q 
as compared to the leading one. 

In this paper we will consider quark 27rDA's, defined as matrix elements of the light-cone 
quark string operator: 



(^(p,k)\—J2q(^,-z}zq(-z)\o) = (p-z)f ^$V^V (2M - 1)(p - 2) , 

N fq J 

= (p • z) I dvH Q (v,£,W 2 ) e iv(p - z) , (6) 



with x = 7 M Xu and 



HQ(v,Z,W 2 ) = V(^,£,W/ 2 ). (7) 



[x,y] denotes a path-ordered exponential [x,y\ = Pexp[igf dt (x — y) M A M (tx + (1 — t)y)] 
which ensures gauge-invariance of the matrix element. In (§) we have introduced both 
types of single-varaiable distributions which can be found in the literature - the distribu- 
tion amplitude $ Q (w, £, W 2 ) and the symmetric form H Q (u, £, W 2 ). As we shall see later, 



calculations in the coordinate-space naturally require to introduce a parametrisation in 
terms of a double-distribution function introduced in [031: 



(2tt(P, k) $(x) [x, -x]xij;{-x)\0) = (P • x) / d[a, £]F(a, /?, ^^e* 01 ^^*^ + 



where 



+(k ■ x) J dpD(p, W 2 )e i0{k - x) + 0(x 2 ) 



1 r0--\ a \) 

d\a,0\ = I da d(3 . (9) 

1 J-(l-|a|) 



The second term in (§) corresponds to the so-called "D-term" introduced in [ 24[ . Com- 
paring matrix element (§) evaluated on the light-cone, x 2 = 0, with parametrisation (|6]) 
one easily obtains a relation between double and 27r distribution (0): 

H Q (u, £, W^ 2 ) = | d[a, /3]P(a, /3, iy 2 )5(a + £/3 - u) + 9(\u\ < sign(£) D(u/£) . (10) 
Note that the matrix element of the pseudovector quark string operator 

M^*)|-jJ-E^)[^-^7b?(-«)|o> = o (11) 

vanishes because of the positive parity of the final pion pair. 

Finally, we recall that all distributions introduced here depend also on a normalization 
scale fi. This dependence is not relevant for our discussion and has been neglected. For a 
sake of simplicity from now on we will also neglect the gauge factors [x, y] and, wherever 
possible, the explicit W - dependence. 

3 Angular momentum analysis of 7*7 — > nn amplitude 

Our goal here is to identify the leading contribution to the amplitude P 1 77r7r . For simplicity 
we will consider the massless quark degrees of freedom only. Taking into account quark 
masses and gluonic degrees of freedom does not change our conclusions qualitatively. A 
convenient coordinate frame is obtained by identifying the null vectors p and z introduced 
in (Q) with a pair n, n* which obeys 

n 2 = n* 2 = 0, n ■ n* = 2, , 

n=(l,0,0,-l),n* = (l,0,0,l). [ ' 

Taking 

z = n, p = Q n* . (13) 

one obtains a collinear collision of 7* and 7, with the final pion pair moving in the 
virtual photon direction. The virtual photon can carry both transverse and longitudinal 
polarizations. For asymptotically large Q 2 the interaction between photons occurs at a 
light-like separation. 

Let us consider first the tree-level contribution. Corresponding space-time configuration 
is depicted on Fig.0. The real photon and the final quark-antiquark pair move along light 



rays defined by = and x + = Xo + X3 = 0, respectively. If both photons 

have the same helicities, the projection of the angular momentum onto the collision axis 
J z = 0. Vector couplings create the qq pair with opposite polarizations, and the angular 
momentum conservation requires that the pair has the orbital angular momentum L z = 0. 
Transition of such a pair into a 2 pions with L z = is described by twist- 2 matrix element 
(!)• 

Longitudinal polarization of the virtual photon leads to a configuration with J z = ±1, see 
Fig.|3]. Due to the angular momentum conservation, the quark-antiquark pair is created 
with the intrinsic orbital angular momentum L z = J z . Subsequently, it evolves into a 
pion pair with L z = ±1. The amplitude must be therefore proportional to k±. On the 
other hand, non-zero quark orbital angular momentum requires that the transition occurs 
through matrix element of a quark twist-3 operator, see the next section, which results in 
a suppression k^/Q as compared with the amplitude corresponding to J z = 0. 
So far we have considered only a point-like interaction of the real photon with the virtual 
one. Another contribution arises when the real photon splits into its quark-antiquark 
component long before its interaction with 7*. Configuration arising if the transition is 
described by twist-2, L z = photon distribution amplitude is shown in Fig.|4]. From di- 
mensional considerations it follows that twist-2 matrix element between a real photon and 
the vacuum is parametrized by a constant carrying dimension one. It is usually defined 
as A — xi^^), where x is the vacuum magnetic susceptibility, and is the quark 

condensate p3[ . Such a contribution would therefore provide a term in y l 77r7r of the order 
of fj/Q- Note, however, that the angular momentum conservation requires that the qq 
pair has to acquire in the hard collision two units of orbital angular momentum. It means 
that twist-2 contribution to this amplitude vanishes to all orders in perturbation theory. 
The situation here is similar to the computation of the amplitude of hard photoproduc- 
tion of transverse vector mesons by longitudinal photon, discussed in [E7|, P8[ . Possible 
higher-twist corrections to this mechanism will be suppressed by at least one power of Q. 
We conclude that the process described in the previous paragraph provides the leading 
contribution to T^ 7 ™'. 

4 Twist-3 amplitude T 1 77r7r in the Wandzura-Wilczek approximation 



In this section we employ technique developed in |26[ to obtain twist-3 contribution to the 
amplitude T^ 7 ™ in the Wandzura-Wilczek approximation i.e., neglecting all quark-gluon 
operators. Our stating object is the amplitude ([|) written in the form: 

T" v = 16i f d 4 xe- ix{q ~ q,) (2n(P, k)\TJ"{x)r{-x)\0) . (14) 



The main contribution to the above integral arises from distances in the vincinity of the 
light cone, x 2 ~ 1/Q 2 - 

Let us first consider a situation where both photons are virtual and their momenta are 
spacelike: 

-q 2 = Q 2 » A QCD , -q' 2 = Q' 2 » A QCD (15) 



To extract the part of the amplitude which corresponds to the case when the one of the two 
photons is longitudial we have to separate contribution which is proportional to transverse 
momentum k±. An analogous problem in DVCS on a pion has been considered recently 
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Figure 1: Leading contribution to the amplitude (Eh in the scaling limit. 



in |T7L P, O, EU1 and [21]. Here, we present a solution obtained using coordinate-space 



technique which is close to the approach of [21 



Let us consider the T-product of two electromagnetic currents. In the vincinity of the 
light cone x 2 ~ the dominant terms are given by diagrams depicted on Fig. (JJ). Their 
contribution can be written as: 



T{J"{x)r{-x)} 



1 1 

16vr 2 ^ 



+ie^\aX\[*P(x)'yal5*P(-x) + ^(-x)7 CT 75^(x)] \ . (16) 



Here we have introduced a shorthand notation: 

\(T — 9pL\9vu + 9na9uX — 9iiv9\v ■ 



(17) 



Our goal now is to determine matrix elements of vector and pseudovector non-local op- 
erators ■0(x)7( T '0(— x) and ^){x)^ a ^^l){—x\ To this end we employ the operator identity 



derived in [19 



-i/;(x)$i/)(—x) + 



(Vx) 

J o 1 da{e^ + e-^} 



2 Jo 



0, 



i J da [e- & 



M (xV) 
V 



(xd) 



^(ax^ip^—ax) + 



{xV) - e^}tt^ k x l7 ^d k ij{ax)^{-ax) + ... 



(Vx) 



where a = 1 — a, <9 M = d/dx^, and ellipses denote contributions of three-point quark-gluon 
operators which are neglected in the WW approximation. Here, T>^ is a derivative with 
respect to translation introduced in [25|: 



d 

V a {i>(tx)r[tx, -tx)ip(-tx)} = — [i)(tx + y)T[tx + y, -tx + y]ip(-tx + y)} 



dy 



(19) 



with a generic Dirac matrix structure V. Hence, when acting on a matrix element between 
vacuum and the final state with a four-momentum P M , — > iP^. Identity allows to 
compute matrix element of ijj(x)^f CT i()(—x) in terms of the matrix element of the symmetric 



operator ip(x)xip(—x), given by (||). Taking matrix elements of both sides of flT8| ) between 
2-pion state and the vacuum one obtains after simple manipulations: 

(2n(P, k)\i;{x)^(-x)\0) = P CT Fi(Px, kx) + k a D x {Px, kx) + 

i[k a (Px) - P a {kx)][F 2 (Px,kx) + D 2 (Px,kx)}. (20) 

Here we have introduced a compact notation 

Fi(Px,kx) = J d[a,(3}F{a,(3)e ia{Px)+il3{kx \ 

F 2 (Px,kx) = - [ d[a,/3]F(a,(3)/3 [ dtt^ + e-^je''^ 1 '^'' 1 ' 1 , 
2 7 Jo 

Dt{kx) = J^df3D(p)e if3{kx) , 

D 2 (Px,kx) = I f 1 d/3D(B) C dtt\S Px) -e-'^le" 3 '^'. (21) 
2 7-i 7o 



Consider now the pseudovector operator ip{x) r j ^^{—x) . Using the identity fl26" 



^(Xb/ZTs^C - ^) = / dt— — ip(tx)^ 5 ip(-tx) - ienuap dttx u T> a ip(tx)^ip(-tx) 

JO OXn JO L 



1 /■* 







dt I dv if)(tx) {tgG^ v {vx) + v r y 5 igG flu (vx)} x u $ip(—tx) , 

(22) 



one can express its matrix elements in the WW approximation through matrix elements 
of the vector operator ip^aip. In this way one obtains 

(2n{P, k)\${x) lal ^{-x)\U) = e aXpS x x P p k 5 [F 3 {Px, kx) + D 3 {Px, kx)} (23) 

where 

F 3 (Px,kx) = - [d[a,0\F(a,P)p t dtt[^ Px ^ - e -^)] e ^W^)+*)] ; 
2 7 Jo 

D 3 (Px,kx) = - j 1 d(3D{(3) t dtt[S Px) + e- il{Px) }e m{kx) . (24) 



2 7-i Xl ' Jo 

Now, combining (P0| ) and (|23|) with ([U|) and (|l~4l) one finally finds 

= ^Jd 4 xe- lx ^'^l [ tZxxPa{F 1 (Px,kx)-F 1 (Px,kx)} + 

+ ix^kl, P a ]{F 2 {Px, kx) - F 3 {Px, kx) + F 2 (Px, kx) + F 3 (Px, kx)} + 
+ ix u x a [k1, P a }{F 2 (Px, kx) + F 3 (Px, kx) + F 2 (Px, kx) - F 3 (Px, kx)} + 

+ ix 2 [k^P u }{F 3 (Px,kx) -F 3 (Px,kx)} + ■••}, (25) 

where for simplicty we have denoted by ellipses all D-type contributions, Fi(Px,kx) = 
Fi(—Px, —kx) and 

[Jtf, P CT ] = MP" - k w ± P" = k»P u - k v P» + 0{W 2 /Q 2 ) . (26) 



The last line follows from the light-cone expansion (H). 

Calculation of the Fourier integrals is now straightforward. Note that at this point the 
advantage of using the double-distribution representation for the matrix element (§) be- 
comes clear, as all dependence on the momentum k^ is in the exponents. Separating the 
linear in k± term one obtains after simple manipulations: 



IF 



£• 



q Q 2 



d[a,p}F(a,p) 



kl{P + uq'ydt:C{a + t(3,w) + 



+ e x {P - uq) v d ( C(a + £/3, 



-w) 



where 



to 



q 2 — q' 2 
q 2 + q' 2 ' 



-1< u < 1, 



and 



C(u, to) 



lo 2 {1-u) 



In 



1 + LOU 



10 



LO 2 1 + u) 



In 



LOU 



l+LO 



(27) 



:28i 



(29) 



Using an identity: 

\ d[a, /3]F(a, p)d i C(a + w) = du d[a, 0\F(a, (3)5{a + £(3- u)C{ 

% j ^duH$(u,OC(u,w), where H$(u,£) = J d[a, 0\F(a, (3)5{a + - 



u, w 



u 



one arrives at 



i if — i 



q Q 2 



dudsH$(u,£) 



k v ± (P + Loq'YC(u, w) + A£(P - ujq) u C(u, -w). 



(30) 



(31) 



Calculation of terms involving D-functions is very similar. We have checked that all the 
D-type contributions can be incorporated together with F-type terms into 



(32) 



with the coefficient function C(u,w) given in ( p9| ) and the 2tt distribution function 
H®(u,£) defined according to fllCf ). 

The coefficient function C(u,w) has a smooth behaviour in the limit w — > 0: 



C(u,u;) = 2m + OH, 

In the limit when photon with momentum g' becomes real, q' 2 —>■ or equivalently w 
one finds: 



(33) 



C(u,w= 1) 



(l-«) 



In 1 



1-ti 



(l+«) 



In 1 - 



l+u 



but 



— Au 

lim C(tt, -w) = ln(l - «;)- — + 0(1 + w) 

U)-vl (1+tt)^ 



(34) 



(35) 



such that the amplitude ( |32|) is singular in this limit. This singularity cancels, however, 
when contribution to the scattering amplitude is evaluated by contraction with polariza- 
tion vectors of photons with helicities A and A', respectively: 



(A, A') =Tr e^A^A') 



(36) 



The second term in (|32| ) can give a non-zero contribution only when a transverse pho- 
ton with momentum q collides with a longitudinal photon with momentum q' . All other 
helicity combinations lead to vanishing contribution due to contraction of transverse po- 
larisation vectors with longitudinal momenta. Note that in the c.m.s. of two photons the 
logitudinal polarisation vector e^O) can be written as ||: 



e„(0) 



-q 



r2 



1/2 



.(qq'f 



q 2 q' 2 



q-q 



(37) 



Evaluating the contraction (|36|) one finds that an additional factor \/l — w arises which 
cancels the logarithmic singularity fl35|). Hence, for the amplitude J 1 ] 77 ™ of a collision of the 
virtual photon with momentum q with a real photon with momentum q', as introduced 
in (|^), one finds: 



or in terms of 2tt distribution amplitude 



In 1 



ln(l 



u 



u 



1 + u) 



ln(n) 
[l-u) 



In 1 



1 + u 



, (38) 



(39) 



Discussion 



Recall that due to the positive charge parity, the pion pair is produced in a configuration 
symmetric with respect to interchange of two pions. Obviously, the production amplitude 
(||) should exhibit the same symmetry. As the part involving 7\ is proportional to k^, 
which is odd under exchange of pions, T x itself must be an odd function of £ as well. 
Note that although H®(u,l;) and $^(m,^) are symmetric functions of £, the presence 
of ^-derivative in expressions (|38| ) and (|39D indeed results in Ti which is antisymmetric 
function of £. 

In order to estimate the magnitude of Ti, we recall the 'minimal' model of the quark 27r 



DA ||, || 



® Q (u, C, W 2 ) = -30u(l - u) (2u - l)—M Q (fi 2 )B(^, W) . 



(40) 



Here M^(n 2 ) is the momentum fraction carried by quarks in a pion at a scale /i 2 . Function 
B(£, W) is related to the so-called Omnes functions JHJ ^(H 7 ) and ^(H 7 ) by 



3C-/3 2 



f (W) P (cos9 c 



P 2 f 2 (W) P 2 (cos9 cm ) 



(41) 



where C — 1 + O(m^). Omnes functions were analyzed in details in [|3T]. fo{W) and 
fiiW) can be obtained from dispersion relations derived in ||: 



/,(W) = exp 



+ Re / ds 

TV J Ami 



s(s-W 2 -iO) 



(42) 



Here <5 ; are tttt scattering phase shifts in the correspondig channels. 

After taking the derivative over £, the fo contribution vanishes ane one finds: 

Tr(W,0 = -f E e ? ^M<V)/ 2 (W/)£. (43) 

As a result, the corresponding contribution to the production amplitude @ has in the 
CM system the angular dependence ~ sin(^ cm ) cos(^ cm )e ±t ^ >cm ~ F^(fl cm , cm ) with / = 2 
and m = ±1. As far as the W dependence is considered, the Omnes function f2{W) is 
small above the threshold, where pions are produced mostly in the S-wave Due to 



/ 2 (1270) resonance, f%iyV) has a peak at W — 1.275 GeV. This range of W is therefore 
most suitable to study Ti experimentally, and the W^-dependence provides an important 
test of the formula ([43]). 

As pointed out in by applying different weighting and averaging procedures one can 
directly compare the QCD predictions for different helicity amplitudes with the 7*7 — > ttti 
data. The analysis of Ti presented here allows to answer the question about the magni- 
tude of the matrix elements of quark-gluon operators, which are not accounted for in the 
WW approximation. By comparing prediction for the helicity-zero amplitude To with the 
data one can constraint a model of 2it quark distribution amplitude. The same model 
can be subsequently used to predict the amplitude Ti according to (p8|) and (p9|) . Any 
discrepancy between such a prediction and the data would hint at a very interesting sce- 
nario where a significant contribution to Ti arises from twist-3 quark-gluon operators. 
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Figure 2: Space-time development of the 7*7 collision. Arrows denote polarizations of photons 
and quarks, respectively. As the final qq pair has no orbital angular momentum, its transition 
into the pion pair is described by a matrix element of twist- 2 light-cone quark string operator. 




Figure 3: Space-time development of the 7*7 collision. Arrows denote polarizations of photons 
and quarks, respectively. The final qq pair carries a unit of orbital angular momentum. Its 
transition into the pion pair is described by a matrix element of a twist-3 operator. 




Figure 4: The real photon splits into a quark-antiquark pair long before a collision with the 
virtual one. Arrows denote polarizations of photons and quarks, respectively. The final qq pair 
has to carry two units of orbital angular momentum and therefore its transition into the pion 
pair must be described by a matrix element of a higher-twist operator. 



